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Abstract — This paper covers recent developments in the 
theory of negative imaginary systems and their application to 
the control of highly resonant flexible structures. The theory 
of negative imaginary systems arose out of a desire to unify a 
number of classical methods for the control of lightly damped 
structures with collocated force actuators and position sensors 
including positive position feedback and integral force feedback. 
The key result is a stability result which shows why these 
methods are guaranteed to yield robust closed loop stability in 
the face of unmodelled spillover dynamics. Related results to be 
presented connect the theory of negative imaginary systems to 
positive real systems theory and a negative imaginary lemma 
has been established which is analogous to the positive real 
lemma. The paper also presents recent controller synthesis 
results based on the theory of negative imaginary systems. 

I. Introduction 

The problem of vibration control for flexible structures 
arises in a variety of applications such as large space struc- 
tures, flexible dynamics in air vehicles, nano-positioning, 
the control of atomic force microscopes and the control of 
optical interfermetric sensing systems; e.g., see [1], [2]. A 
particular problem in vibration control for flexible structures 
is that unmodelled (spillover) dynamics can severely degrade 
the control system performance or lead to instability. Also, 
uncertainties in mode frequencies and damping can lead 
to similar problems of poor control system performance or 
instability; e.g., see [3], [4]. 

An approach to the control of flexible structures which 
addresses these robustness issues is positive-position feed- 
back, which involves the use of collocated force actuators 
and position sensors [5], [6]. Roughly speaking a positive 
position feedback controller consists of a transfer function 
matrix which has the same form as a finite mode model of 
a flexible structure. This approach was originally proposed 
as an alternative to passivity based velocity feedback ap- 
proaches in the control of large space structures. The orig- 
inal papers on positive position feedback included stability 
proofs based on Routh-Hurwitz and Nyquist techniques but 
the robustness properties were never fully quantified. The 
research presented in this paper aims to develop a general 
systems theory framework for positive position control and 
related control methods for flexible structures. An important 
special case of positive-position feedback involves the use 
of piezoelectric actuators and sensors. A positive-position 
feedback controller can be designed to increase the damping 
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of the modes of a flexible structure. Furthermore, this con- 
troller is robust against uncertainty in the modal frequencies 
as well as unmodelled plant dynamics. We consider the 
theory of negative-imaginary systems to reveal the robustness 
properties of MIMO positive-position feedback controllers 
and related types of controllers for flexible structures; e.g., 
see [1], [7]-[15]. 

The negative-imaginary property of linear systems can be 
extended to nonlinear systems through the notion of coun- 
terclockwise input-output dynamics [16]— [18]. It is shown 
in [19] for the single-input, single-output (SISO) linear case 
that the results of [16], [17] guarantee the stability of a 
positive-position feedback control system in the presence 
of unmodelled dynamics and parameter uncertainties that 
maintain the negative-imaginary property of the plant. 

A. Flexible Structure Modeling 

In modeling an undamped flexible structure with collo- 
cated force actuators and position sensors, modal analysis 
can be applied to the relevant partial differential equation 
leading to the transfer function matrix 



i 



S 2 + KiS + W, 



where, for all i, m > 0, u>i > 0, and ipi is an m x 1 vector. 
Therefore, the He rmitian- imaginary part 

$SK[P(JU)] = ~](P(JU)-P*(ju)) 

of the frequency response matrix P(juj) satisfies 



9fH[P0w)] = -wX] 



+ UJ 2 K 2 



for all lu > 0. That is, the frequency response matrix 
has a negative-semidefinite Hermitian-imaginary part for all 
uj > 0. We thus refer to the transfer function matrix P(s) as 
negative imaginary. Any flexible structure with collocated 
force actuators and position sensors will have a negative 
imaginary transfer function matrix. 

II. Negative Imaginary Transfer Function 
Matrices 

In this section, we present the formal definition of a 
negative imaginary transfer function matrix. We also present 
a number of extensions to this notion, which have appeared 
in the literature. 



Definition 1 ( [20]-[23]): A square real-rational proper 
transfer function matrix R(s) is termed negative imaginary 
(NI) if 

1) R(s) has no poles at the origin and in 5ft [,s] > 0; 

2) j[R(jw)—R*(jw)] > for all lo G (0, oo) except values 
of u> where jtu is a pole of R{s)\ 

3) If jwo, loq G (0, oo), is a pole of R(s), it is at most a 
simple pole, and the residue matrix Kq = lim s _> ja , (s — 
juJo)jR(s) is positive semidefinite Hermitian. 

In the SISO case, a transfer function is negative imaginary if 
and only if it has no poles in open right half plane (ORHP) 
or at the origin and its phase is in the interval [— ir, 0] at 
all frequencies. Consequently, the Nyquist plot of a SISO 
negative-imaginary transfer function lies below the real axis; 
see Figure Q] 
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Fig. 1 . Nyquist plot of a negative-imaginary system. 

Remark 1: Note that this definition has been recently 
extended to allow for a simple pole at the origin; [24], [25]. 
However, in this case, this leads to a slightly different form 
of the stability condition. Also, in [26], a lossless version of 
this definition is given along with a corresponding negative 
imaginary lemma. 

The following lemma from [21], [23] describes the rela- 
tionship between negative imaginary transfer function matri- 
ces and positive real transfer function matrices. 

Lemma 1 ( [21], [23]): Given a square real-rational 
strictly proper transfer function matrix R(s). Then R(s) is 
negative imaginary if and only if 

1) R(s) has no poles at the origin; 

2) F(s) = sR(s) is positive real. 

The following definition defines the strict version of the 
negative imaginary property. 

Definition 2 (See [20], [21].): A square real-rational 
proper transfer function matrix R{s) is termed strictly 
negative imaginary (SNI) if 

1) R(s) has no poles in 5ft[s] > 0; 

2) j[R(ju>) - R*(juj)] > for lo G (0, oo). 

Remark 2: In the paper [27], the strict negative imaginary 
property is extended to a notion of strongly strict negative 



imaginary (SSNI) systems by imposing the additional con- 
ditions that 

lim ju(R(jw) - R(ju)*) > 

and 

]imj-(R(joj)-R(ju>)*)>0. 

VO LO 

This extension is useful in characterizing the strict negative 
imaginary property using LMIs and in controller synthesis. 

The following lemmas state some useful properties of 
negative imaginary transfer function matrices. 

Lemma 2 ( [21], [23]): A square real-rational proper 
transfer function matrix R(s) is negative imaginary if and 
only if R(oo) = R T (oo) and R(s) = R(s) - R(oo) is 
negative imaginary. 

Lemma 3 ( [15]): If the m x m transfer function matrix 
R(s) is NI and the m x m transfer function matrix R s (s) is 
NI, then the transfer function matrix 

R T (s) = R(s) + R s (s) (1) 

is NI. Also, if in addition R s (s) is SNI, then R T {s) is SNI. 

Lemma 4 ( [15]): Consider the NI transfer function ma- 
trices R{s) and R s (s), and suppose that the positive- 
feedback interconnection of this transfer functions is inter- 
nally stable. Then the corresponding 2m x 2m closed-loop 
transfer function matrix 

r i?( s )(/-i? s ( s ) J R( s ))- 1 

1 ; L R s (s) (I - R(s)R s (s))- 1 R(s) 

i?( s )(/-i? s ( s ) J R( s ))- 1 i? s ( s ) 1 

R s (s)(I-R(s)R s (s)y 1 \ K) 

is NI. Furthermore, if in addition, either R(s) or R s (s) is 
SNI, then T(s) is SNI. 

We can also can also extend the above definitions of NI 
systems to the finite bandwidth case. 

Definition 3 ( [28]): A square real-rational proper trans- 
fer function matrix R(s) is said to be finite frequency 
negative imaginary (FFNI) with bandwidth Co if it satisfies 
the following conditions: 

1) R(s) has no poles at the origin and in the open right- 
half of the complex plane; 

2) j[R(jio) - R*(jLo)} > for all lo G fl, where 

Q = {lo G K : < lo < lo, jlo is not a pole of R(s)}; 

3) Every pole of R(s) on jCl, if any, is simple and the cor- 
responding residue matrix of jR(s) is positive semidef- 
inite Hermitian, where O = {ojgIR:0<cj< Co}; 

4) R(oo) = R T (oo). 

Based on the above definition, we have the following 
properties. 

Lemma 5 ( [28]): If R(s) is an FFNI transfer function 
matrix with bandwidth Co, then the following properties hold: 

1) j[R(jLo) - R*(jLo)} < for all lo G fii, where 

fli = {lo G K : —Co < lo < 0, jlo is not a pole of R(s)}. 

2) Every pole of R(s) in jfli, if any, is simple and the cor- 
responding residue matrix of jR(s) is negative semidef- 
inite Hermitian, where J7i = {a; G K : — Co < lo < 0}. 



Note that the stability result given in [29], uses a notion 
similar to the FFNI property defined above along with a small 
gain property to establish robust stability. 

III. Some Negative Imaginary Lemmas 

The following Theorem provides a state space condition 
for a system to be negative imaginary. 

Theorem 1 (Negative Imaginary Lemma, [20]-[23].): 
Let (A, B, C, D) be a minimal state-space realization of an 
to x to real-rational proper transfer function matrix R{s), 
where A G M" xn , B G R" x "\ C G R mx, \ D G R raXm . 
Then R{s) is negative imaginary if and only if 

1) det(A) + 0, D = D T , and 

2) there exists a matrix Y = Y T > 0, Y G R nxn , such 
that 

AY + YA T < 0, and B + AYC T = 0. 
The following alternative form of the negative imaginary 
lemma also holds in the case that the system has a simple 
pole at the origin. 

Theorem 2 (Modified NI Lemma, [22]-[25].): Let 
(A, B, C, D) be a minimal realization of R(s) G K mxm , 
where A G M' iX '\ B G R" x "\ C G R mx ", P G R TOXm , 
to < n. Then R(s) is negative-imaginary if and only if 



D 

W G K mx ' 
is satisfied: 



D T and there exist matrices P = P T > 0, P G R nXT \ 
R mxm , and L G R mx ™ such that the following LMI 



PA 4 
P T P 



A T P 
-CL4 - 

—L T L 
—W T L 



PB - A T C T 
(CB + B T C T ) 

-L T W 
~W T W 



< 0. 



(3) 



Remark 3: Note that a version of this negative imaginary 
lemma for descriptor systems is presented in [30]. 

The following theorems provide corresponding conditions 
for strict negative imaginary properties. 

Theorem 3 (Weak SNI Lemma, [22], [23].): Let 
(A,B,C,D) be a minimal state-space realization of 
an to x to real-rational proper transfer function matrix R(s), 
where A e R" xn , B e E nxm , C G R mxn , D e W nxm . 
Then R(s) is strictly negative imaginary if and only if 

1) A is Hurwitz, D = D T , rank(P) = rank(C) = m; 

2) There exists a matrix Y = Y T > 0, Y G R" xn , such 
that 

AY + YA T < 0, and B + AYC T = 0; 

3) The transfer function matrix M(s) ~ 



.4 



LY~ L A~ 



B 







has full column rank at s = joj 

for any u £ (0, oo). Here L T L = -AY - YA T . That 
is, rank(Af(jo;)) = to for any uj G (0, oo). 
The following theorem provides a purely LMI approach 
to test the strong SNI property. 

Theorem 4: [Strong SNI Lemma, [27], [31].] Given a 
square transfer function matrix R(s) G TZ mxm with a 
state-space realization (A, B, C, D), where A G IR™ xn , 
B G M" xm , C G R mx " and P G M mXm . Suppose 
R{s) + R{— s) T has normal rank m and (C, A) is observable. 



Then, A is Hurwitz and R(s) is Strongly Strict Negative 
Imaginary if and only if D = D T and there exists a matrix 
Y = Y T > such that 



AY + YA T < and B = —AYC T . 



(4) 



The following theorem from [28] gives a finite bandwidth 
version of the Negative Imaginary Lemma. 

Theorem 5 (FFNI Lemma [28].): 
Consider a real-rational proper 

transfer function matrix R(s) with a minimal state- 
space realization (A, B, C, D). Suppose all poles of R(s) 
are in the closed left-half of the complex plane, and the 
poles on the imaginary axis, if any, are simple. Let a 
positive scalar Q be given. Also, suppose that if A has 
eigenvalues joji (i G {!,..., q}) such that < < uj, 



the residue of (si — A) at s 



jwi is given by 



lim.. 



(s — juji){sl — A) . Then the following 



statements are equivalent: 

1) The transfer function matrix R(s) is FFNI with band- 
width UJ. 

2) det(A) ^ 0, D = D T , and the transfer function 
matrix F(s) with a minimal state-space realization 
(A, B, CA, CB) is FFPR with bandwidth Q. 

3) det(A) ^ 0, D = D T , and CA^iB = (CA^B)* > 
for all i G {1, . . . , q} if A has any eigenvalues on jfi 
where f2 = {a;G]R:0<w< ui}. Also, there exist real 
symmetric matrices P = P T and Q = Q T > such 
that 

PA + A T P - A T QA + u~! 2 Q < 0, (5) 
C + B T A- T P = 0, (6) 
QA~*B = 0. (7) 

4) det(A) / 0, P = P T , and CA^iB = (CA^B)* > 
for all i G {1, . . . , q} if A has any eigenvalues on jfi 
where f2 = {a;G]R:0<w< a)}. Also, there exist real 
symmetric matrices Y = Y T and X = X T > such 
that 

AY + YA T - AXA T + uj 2 X < 0, (8) 
B + AYC T = 0, (9) 
CX = 0. (10) 

In addition to the above negative imaginary lemmas, which 
involve LMI conditions to test if a given system is negative 
imaginary or strictly negative imaginary, the following result 
uses an eigenvalue test of a corresponding Hamiltonian 
matrix to test if a given system is NI. 

Theorem 6 ( [32]): Consider a transfer function matrix 
R(s) with a minimal realization (A, B, C, D). Also, suppose 
A is a Hurwitz matrix, Q = CB+B T C T > 0, and D = D T . 
Then, R(s) is NI if and only if the following conditions are 
satisfied: 

1. The Hamiltonian matrix, 

PQ-!P T 

A T _ A T C Tq-1 B T 



N 



A + BQ 
-A T C T Q 



~ X CA 
~ X CA 



has no pure imaginary eigenvalues with odd multiplicity. 



2. jif^ii±!i±ll (R^jTjj) — R(jr]i)*) has only positive real 
eigenvalues for all i G {l,fc — 1}, where u>i G 
il, and il = {lji,u>2, ■ ■ ■ ,Wk} is me set of fre- 
quencies listed in strictly increasing order, such that 
det[F(jwi) +F(juji)*] = 0. Here, F(s) = sR{s) 
where R(s) = R(s) - D. 
Remark 4: This approach has been extended in the paper 
[33] to consider the problem of perturbing a non-NI system 
model in order to obtain an NI system model. 

Also, the paper [32] gives the following result for testing 
the NI property of SISO systems. 

Assumption 1: Let A G W lxn ,B G R" xl ,C G 
M. lxn ,D G M, and suppose the transfer function G(.s) = 
C(sl — A)~ 1 B + D has all of its poles and zeros in the 
closed left half of the complex plane excluding poles at the 
origin. Also, any pole on the imaginary axis is assumed to 
be simple. 

Theorem 7 ( [32]): The SISO transfer function G(s) = 
C(sI-A)~ 1 B+D with CB > 0, and satisfying Assumption 
Q]is NI if and only if the following conditions are satisfied: 

1. The matrix (I - T ^BC)A 2 has no eigenvalues of odd 
(algebraic) multiplicity on the open negative real axis. 

2. All residues of F(s) — s (G(s) — D) corresponding to 
poles on the imaginary axis are positive. 

IV. Stability of Interconnections of Negative 
Imaginary Systems 

The following theorem, which is the main stability result in 
the theory of negative imaginary systems, gives a necessary 
and sufficient condition for the stability of the positive feed- 
back interconnection (see Figure |2]i between an NI system 
and an SNI system. 
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Fig. 2. A positive feedback interconnection. 



Theorem 8 ( [20]-[23]): Given a negative imaginary 
transfer function matrix R(s) and a strictly negative imag- 
inary transfer function matrix R s (s) that also satisfy 
i?(oo)i? s (oo) = and R s (oo) > 0. Then the positive 
feedback interconnection [R(s), R s (s)] is internally stable 
if and only if A max (i?(0)i? s (0)) < 1. 

Remark 5: Note that in the papers [34], [35], the robust- 
ness properties obtained by this theorem are compared with 
the robustness properties which are obtained from the [i 
analysis approach. Also, in the papers [36], [37], this stability 
robustness result is extended to allow for additional norm 
bounded uncertain parameters. 



The above theorem can also be extended to give a suffi- 
cient condition for closed loop stability in the case in which 
the NI system has a simple pole at the origin. 

Theorem 9 ( [24], [25]): Consider a strictly proper neg- 
ative imaginary transfer function matrix R(s) with minimal 
state space realization (A,B,C, 0), which may contain a 
simple pole at the origin. Also, consider a strictly negative 
imaginary transfer function matrix R s (s). Then the closed- 
loop positive feedback interconnection between R(s) and 
R s (s) is internally stable if R s (s) < and the matrix 
A + BR s (0)C is not singular. 

V. Negative-Imaginary Feedback Controllers 

In this section, we present a number of different classes 
of commonly used negative imaginary controllers; see also 
[15]. 

A. Positive-Position Feedback Controllers 

In the SISO case, a positive-position feedback controller 
(e.g., see [5], [6], [12], [38]-[44]) is a controller of the form 

2—1 l 

where u; 4 > 0, Q > 0, and ki > for i = 1, 2, . . . , M. Such 
controllers are always SNI. 

B. Resonant Controllers 

Resonant controllers are proper but not strictly proper 
SISO SNI controllers of the form 

M -h 2 

°M-£ ?T^Tsy ,12) 

where > 0, Q > 0, and ki > for i = 1, 2, . . . , M; e.g, 
see [8], [45]. The controller ([T2l can be implemented as the 
positive-position feedback controller (ITT1 using an accelera- 
tion sensor rather than a position sensor. Alternatively, ([T2l 
can be implemented as the positive-real feedback controller 

2—1 1 

where uj, > 0, Q > 0, and fc 4 > for i = 1,2, ...,M, 
using a velocity sensor rather than a position sensor. Such 
controllers are always SNI. 

C. Integral Resonant Controllers 

MIMO transfer function matrices of the form 

c(s) = [si + r^r 

are SNI. Here, T is a positive-definite matrix and 4> is a 
positive-definite matrix. The use of a controller of this form 
when applied to a flexible structure with force actuators and 
position sensors is referred to as integral resonant control, 
or integral force control [1], [9], [10], [46]-[51]. 



VI. State-Feedback Controller Synthesis 



VII. Conclusions 



Consider the feedback control system in Figure [3] in the 
case that full state feedback is available. In this case, the 
negative imaginary lemma can be used to synthesize a state- 
feedback control law such that the resulting closed-loop 
system is NI. Indeed, suppose the uncertain system shown 
in Figure [3] is described by the state equations 



x = Ax + B\w + B 2 u, 
z = C\x, 
w = A(s)z, 



(13) 
(14) 
(15) 



where the uncertainty transfer function matrix A(s) is as- 
sumed to be SNI with |A max (A(0))| < 1 and A (00) > 0. 
Applying the state-feedback control law u = Kx yields the 
closed-loop uncertain system 



x = (A + B 2 K)x + Biw, 
z = C±x, 
w = A(s)z. 



(16) 
(17) 
(18) 



The corresponding nominal closed-loop transfer function 
matrix is 



G d (s) = Ci(sl — A — B 2 K)- x Bx. 



(19) 



k>J Ll L I ILL I Clt^L'J Ll'Ktp 

G c i(s) 



Fig. 3. A feedback control system. 

Theorem 10 ( [15], [52]): Consider the uncertain system 
( Tol l. (HP , (fl"5"l i and suppose there exist matrices Y > 0, M, 
and a scalar e > such that 



AY + YA T + B 2 M + M T 
Bj + dYA T + CiM 



+ el 



<0, 



AYCJ 




B 2 MCl 



dYCl - I < 0, 
Y > 0. 



(20) 
(21) 
(22) 



Here the parameter e > is chosen to be sufficiently small. 
Then the state-feedback control law u = MY _1 x is robustly 
stabilizing for the uncertain system (Qj), (fl4l . dl5l . 

Remark 6: Note that in the papers [53], [54], an alter- 
native approach to controller synthesis is considered by 
connecting the negative imaginary synthesis problem to the 
H°° control problem. 



We have surveyed some recent results in the area of 
negative imaginary systems theory. This results characterize 
the class of negative imaginary systems and give methods for 
testing if a given system is negative imaginary. These results 
also show that the class of NI systems yields a robust stability 
analysis result, which broadly speaking can be captured by 
saying that if one system is negative imaginary and the 
other system is strictly negative imaginary, then a necessary 
and sufficient condition for internal stability of the positive- 
feedback interconnection of the two systems is that the dc 
loop gain is less than unity. This result provides a framework 
for the analysis of robust stability of lightly damped flexible 
structures with unmodelled dynamics. 
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